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ABSTRACT 

High-magnetic-field radio pulsars are important transition objects for understanding the connection 
between magnetars and conventional radio pulsars. We present a detailed study of the young radio 
pulsar J1119— 6127, which has a characteristic age of 1900 yr and a spin-down-inferred magnetic field 
of 4.1 X IQ-'^^ G, and its associated supernova remnant G292.2— 0.5, using deep XMM-Newton and 
Chandra X-ray Observatory exposures of over 120 ks from each telescope. The pulsar emission shows 
strong modulation below 2.5 keV, with a single-peaked profile and a large pulsed fraction of 0.48±0.12. 
Employing a magnetic, partially ionized hydrogen atmosphere model, we find that the observed pulse 
profile can be produced by a single hot spot of temperature 0.13 keV covering about one third of the 
stellar surface, and we place an upper limit of 0.08 keV for an antipodal hot spot with the same area. 
The nonuniform surface temperature distribution could be the result of anisotropic heat conduction 
under a strong magnetic field, and a single-peaked profile seems common among high-i? radio pulsars. 
For the associated remnant G292.2— 0.5, its large diameter could be attributed to fast expansion in a 
low-density wind cavity, likely formed by a Wolf-Rayet progenitor, similar to two other high-i? radio 
pulsars. 

Subject headings: ISM: individual objects (G292.2— 0.5) — ISM: supernova remnants — pulsars: in- 
dividual (PSR J1119-6127) — stars: neutron — X-rays: ISM 



1. INTRODUCTION 

Over the past two decades, our understanding of 
neutron stars has been revolutionized due to discover- 
ies of several new classes of objects (see iKaspil I2010L 
for a review). An extreme class is magnetars, which 
typically have high spin-down-inferred magnetic fieldsQ 
of 1 0^^-10^^ G and show violent bursting activities 
(see iRea fc"Esposit o 2011; Mcrcghctti 2008; iNg et al.l 
1201 It iScholz et al.l l20l2[ ). It is generally believed 
that their X-ray emission i s powered by the decay of 
ultra-strong magnetic fields (iDuncan fc Thompson|[T992l: 
[Thompson fc DuncanI [19951 I1996D . " However, the ori- 
gin of magnetars and their relation to the more con- 
ventional rotation-powered pulsars remains a puzzle. 
The distinction between these two classes of objects 
became increasingly blurred thanks to the recent dis- 
coveries of a relatively low S-field magnetar of -B < 
7.5 X 10^2 G HRea et al.il2010D. magnetar-like bursts from 
a rota t ion- powered pulsar ("Gavriil "erani200llNg et al.l 
[2008at Eimar fc Safi-Harb 200 17 and radio emission 
from magnetars (jCamilo et'aLl 120061 120071 : ILevin et al.l 
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The dipole S-field can be estimated by B = 3.2 X 
10l9(pp)i/2 G, where P is the spin period in second and P is 
the spin-down rate. 



[20T0I) . These provide some hints that magnetars and 
radio pulsars could be drawn from the same population, 
and the former could represe nt the high-field tail of a sin- 
gle birth g -ficld distribution (jKaspi fc McLaughlinll2005l: 
iPerna fc P ons 2011; Ho 2012). One way to test this "uni- 
fication" picture is via observations of high-magnetic- 
field radio pulsars^ which are a critical group of rotation- 
powered pulsars that have similar spin properties as mag- 
net a.rs, imp lying magnetar-like dipole B-field strengths 
(see iNg fc K aspi 20111 for a review). These pulsars could 
represent transition objects that share similar observa- 
tional properties with magnetars and radio pulsars, pro- 
viding the key to understanding the connection between 
the two classes. 

PSR J1119— 6127 is a young high-i? radio pulsar 
discovered in the Parkes Multibeam Pulsar Survey 
([Camilo et al. 2000). It has a spin period P — 408 ms 
and a large period derivative P = 4.0 x lO"^'^, which to- 
gether imply a surface dipole i?- field of 4.1 x 10^^ G. It is 
one o f the very few pul s ars w ith a measured braking in- 
dex^ [WtevrSeeLaD (llOp) reported n = 2.684±0.002 
using over 12 years of radio timing data, which give a 
characteristic age of P/[(n—l)f'] = 1.9 kyr. The true age 
could even be younger if the pulsar were born with a spin 
period close to the present-day value and n is constant. 
The pulsar is associated with a 17'-diamctcr radio super- 
nova remnant (SNR) shell G292.2-0.5 (Crawford_eialJ 
120011 ) . Based on Hi abs orption measurements of the 
SNR, iCaswell et al.l ()2004[ ) reported a source distance of 

^ Throughout this paper, we refer to rotation-powered pulsars as 
radio pulsars, even tho ugh their radio beams may miss the Earth, 
e.g. PSR J1846-0258 H Archibald et al.|[200l) . 

^ The braking index is defined as n = i>i} / , where u is the spin 
frequency and u and i> are its first and second time derivatives, 
respectively. 
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TABLE 1 

X-Ray Observations of PSR J1119-6127 Used in This Study 



ObsID 



Date 



Instrument / Mode Net Exp.'' (ks) 



XMM-Newton 
0150790101 2003 Jun 26 PN / Large Window 



0672790101 2011 Jun 14 



0672790201 2011 Jun 30 



41.8 

MOSl / Full Frame 44.3 

MOS2 / Full Frame 44.8 

PN / Large Window 27.7 

MOSl / Full Frame 34.1 

MOS2 / Full Frame 34.3 

PN / Large Window 27.9 

MOSl / Full Frame 33.8 

MOS2 / Full Frame 33.7 



Chandra 


2833 


2002 Mar 31 


ACIS-S / TE 


56.8 


4676 


2004 Oct 31 


ACIS-S / TE 


60.5 


6153 


2004 Nov 02 


ACIS-S / TE 


18.9 



After removing periods of high background. 

8.4±0.4kpc. 

In the X-ray band, the pulsar and SNR were de- 
tected with ASCA, ROSAT, t he Chandra X-ray Ob 
servatory, and XMM-Newton (Pivovaroff et al. 20011 
Gonzalez fc Safi-Harb 2003, 2005; Gonzalez et al. 200i 
Safi-Harb fc Kumaiil2008 il. The pulsar e mission consists 



of thermal and non-thermal com ponents ([Gonzalez et alj 
120051 : iSafi-Harb fc Kumadl200l ) and shows strong pu l- 
sations in the soft X-ray band ([Gonzalez et al.l [20051 ) . 
There is a faint, jet-like pulsar wind nebula (PWN) 
extending 7" from the pulsar ([Gonzalez fc Safi-Harbl 
[2003; Sa fi-Harb fc Kumaiil2008D. but n o radio counter- 
part has been found ([Crawford et all 120011) . The in- 
terior of SNR G292. 2-0.5 shows faint, diffuse X-ray 
emis sion with spectral lines that indicate a thermal ori- 
gin ([Gonzalez fc Safi-Harbl [2"005| ) . Recently, the Fermi 
Large Area Telescope detected 7-ray pulsations from 
PSR J1119-6127, making it the highe st B-field 7-ray 
pulsar yet known ([Parent et al.l 120 lit) . Here we re- 
port on a detailed study of PSR J1119-6127 and SNR 
G292.2— 0.5 using new XMM-Newton observations to- 
gether with archival XMM-Newton and Chandra data. 

2. OBSERVATIONS AND DATA REDUCTION 

New observations of PSR J1119— 6127 were carried out 
on 2011 June 14 and 30 with XMM-Newton. The PN 
camera was operated in large window mode, which has 
a time resolution of 48 ms, while the MOS cameras were 
in full-frame mode, with 2.6 s time resolution, both with 
the medium filter. We also reprocessed archival XMM- 
Newton and Chandra data taken in 2002-2004. The 
Chandra observations were all made with the ACIS-S 
detector, which has 3.2 s frame time. Table [T] lists the 
observation parameters of all data sets used in this study. 

XMM-Newton data reduction was performed using 
SAS 11.0 with the latest calibration files. We pro- 
cessed the Observation Data Files with the tools 
emproc/epproc, then employed espfilt to remove pe- 
riods of high background. The resulting exposure for 
the PN and MOS cameras are 97 and 112 ks, respec- 
tively. Throughout the analysis below, we used events 
having PATTERN < 12, unless specified otherwise. Chan- 
dra data reduction was carried out with CI AO 4.4 and 
CALDB 4.4.7. We applied the script chandra_repro to 
generate new level=2 event files with the latest time- 
dependent gain, charge transfer inefficiency correction 



and sub-pixel adjustments. We obtained a net exposure 
of 136 ks after rejecting periods with background flares. 

3. ANALYSIS AND RESULTS 

3.1. Imaging 

We generated individual exposure-corrected images for 
each PN and MOS observation with the tasks evselect 
and eexpmap, then combined them using emosaic. Fig- 
ure [1] shows a three-color XMM- Newton image overlaid 
with 1.4 GHz radio contours from lCrawford et al.l ([20011 ) 
and a broadband intensity image in the 0.3-8 keV energy 
range. We have also generated images with the Chandra 
data, but they are not shown here due to the poorer sen- 
sitivity and sky coverage. As shown in the figure, PSR 
J1119-6127 and a partial sheU of SNR G292. 2-0.5 are 
detected. The SNR emission is clumpy and it generally 
traces the brightest part of the 17' diameter radio shell, 
except for the diffuse emission near the pulsar, which 
shows no radio counterpart. The western half of the shell 
is brighter and redder than the eastern half. Note that 
the southwest rim of the shell falls outside the PN camera 
field of view, hence, the sensitivity is largely reduced. 

3.2. Timing 

The timing analysis was carried out with the PN data 
only, since only they have sufficiently high-time resolu- 
tion. We extracted source events from a circular aperture 
of 18" radius centered on the pulsar. This radius was 
determined by the tool eregionanalyse to yi eld the op- 
timal signal-to- noise ratio. We followed Gonzalez et alJ 
([20051 ) to correct for a timing anomaly in the 2003 data; 
we did not find any other such anomaly in the other ob- 
servations. We divided the data into low-energy (0.5- 
2.5 keV) and high-energy (2.5-8 keV) bands, and ob- 
tained 889 ± 30 and 414 ± 20 total counts, respectively, 
with ~ 110 and 120 background counts, respectively, as 
estimated from nearby regions. The event arrival times 
were first corrected to the solar system bary center, then 
folded according to the radio ephem erides. We used 
the p ublished pulsar ephemeris from IWeltevrede et al] 
([20111 ) for the 2003 data, while the ephemeris in 2011 
was obtained from contemporaneous radio observations 
with the Parkes ra dio telescope, as part o f the timing 
program for Fermi ([Weltevrede et al.ll2010D . 

First, we verified that the pulse profiles in 2003 and 
2011 were consistent. They were subsequently co-added 
in the analysis to boost the signal-to-noise ratio. The 
resulting profile is shown in Figure [21 The pulsar emis- 
sion below 2.5 keV is highly modulated. It exhibits 
a single-peaked profile with good alignment with the 
radio pulse. A sinusoid provides an adequate fit to 
the profile and it gives an rms pulsed fraction (PF) of 
0.3 8 lb 0.03. To compare with previous work, we fol- 
low I Gonzalez et al.l ([2005D and calculate the "max-min 

PF" using PF=(F„,ax-^;„in)/(Fmax + F,ni„), whcrC 

and i^niin are the maximum and minimum background- 
subtracted counts in the binned profile, respectively. In 
this way, we obtain max-min PF of 0.48 ± 0.12 in the 
0.5-2.5 keV range. This is somewhat low e r than the value 
0.74 ± 0.14 reported bv iGonzalez et al.l ([20051 ). and the 
discrepancy could be attributed to Poisson fluctuations 
or to difference in phase binning. On the other hand, 
no pulsations are detected above 2. 5 keV. T his is con- 
sistent with the findings of .Gonzalez et al.l ([2005i ) and 
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Fig. 1.— Three-color image of PSR J1119-6127 and G292. 2-0.5 in the 0.3-1.5 keV (red), 1.5-3 keV (green), and 3-8 koV (blue) bands, 
made with all XMM-Newton PN and MOS data listed in Table [T] with a total exposure of 100 ks. The image is exposure-corrected, 
smoothed to 13" resolution, and overlaid with 1.4-GHz radio contours from [Crawford et al.l 1120011 ). Inset: intensity map in the 0.3-8 keV 
band, indicating the source and background extraction regions for the spectral analysis. 

3.3.1. Phase-averaged Pulsar Spectroscopy 



Spectra of PSR J1119-6127 were extracted from 18"- 
radius apertures from the XMM-Newton and Chandra 
data. Background spectra were extracted from nearby 
regions free from SNR emission shown in Figure [TJ We 
have also tried using annular background regions cen- 
tered on the pulsar, and the results are very similar. 
There are 1062 ± 37, 500 ± 25, 442 ± 24, and 908 ± 35 net 
pulsar counts in 0.5-8 keV from the PN, MOSl, M0S2, 
and ACIS detectors, respectively. A comparison of differ- 
ent spectra shows no variability between epochs. There- 
fore, we performed a joint fit to all the data. Given 
the low number of source counts, we did not attempt 
to account for the cross-calibration uncertainty between 



wi th the RXTE result s in the 2-60 keV range reported 
bv [Parent et al.l (|2011[ ). We place put an upper limit of 
PF < 10% in the 2.5-8 keV range. 

3.3. Spectroscopy 

We carried out the spectral analysis with XSPECP^ 
12.7.1. Spectra were grouped to at least 20 counts per 
bin, and only single and double events (i.e. PATTERN < 
4) were used for the PN data to ensure the best spec- 
tral calibration. We restricted the fits in t he 0.5-8 keV. 
A pho toelectric absorption model tbabs bv I Wilms et al.l 
()2000[ ) was employed throughout our study, and we 
adopted the abundances given by the same authors. 

|http: //heasarc .gsf c.nasa. gov/xanadu/xspec/| 
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Fig. 2.— X-ray pulse profiles of PSR J1119-6127 from XMM- 
Newton FN observations in 0.5-2.5 keV (red) and 2.5-8 keV (blue). 
The solid line shows the model lightcurve of the soft X-ray emis- 
sion from a single hot spot with a magnetized atmosphere (see Sec- 
tion UTTTJ . Background levels are indicated by the dashed lines. 
Phase corresponds to the radio pulse at 1.4 GHz shown in the 
lower panel. The radio flux density is in an arbitrary unit. 

instruments (see iTsuiimoto et al.ll201ll ). 

To account for the PWN contamination, we employed 
a PWN component in the spectral model, which was ex- 
tracted from the Chandra data using an annular aperture 
of 2.5"-18". We fitted the PWN spectrum using an ab- 
sorbed power-law (PL) model with the column density 
(TVh) fixed at the pulsar value. The PWN component 
was then held fixed during the fit of the pulsar spec- 
trum. Once a new value of A^h was determined, we re-fit 
the PWN spectrum and the whole process was carried 
out iteratively until consistent spectral parameters were 
found. We first tried simple models, including an ab- 
sorbed blackbody (BB) and an absorbed PL for the pul- 
sar spectrum. However, they fail to describe the data, 
and we found that an absorbed BB-f PL model provides 
a much better fit {xt — 0.97). This gives a high effective 
temperature of kT — 0.21 ± 0.04 keV with a emission ra- 
dius of Sl^Jkm (for a source distance of 8.4 kpc), and a 
PL photon index of 2.1 ± 0.8. The BB component has an 
absorbed flux of 1.2±0.3 x lO-^** ergcm"^ in the 0.5- 
8 keV energy range, and the non-thermal components 
have absorbed fluxes of 2.3 ± 0.4 x lO-^^ergcm"^ s"^ 
and 2.8 x 10^^'' ergcm"^ s"^ from the pulsar and the 
PWN, respectively. The best-fit parameters are listed in 
Table [2] and the PN spectrum is shown in Figure [3l^a). 

Following previous studies, we also tried fitting the 
pulsar spectra with a ne utron star atmosphere model 
(NSA; iZavlin et all 119961 ) plus a PL model. The for- 
mer describes a fully ionized hydrogen atmosphere in 
a i?-field of 10^'^ G. During the fits, the stellar mass 
was fixed at 1.4Mq and we assumed uniform thermal 
emission over the entire surface of 13 km radius. This 
yields a fit that is equally as good as the BB-I-PL model. 



TABLE 2 

Best-fit Blackbody Plus Power-law (BB+PL) and 
Neutron Star Atmosphere Plus Power-law 
(NSA+PL) Models to the Phase-averaged 
Spectrum of PSR ,]1119-6127 



Parameter 




BB+PL 


NSA'^+PL 


ivH cm J 




r, 9 + 0.5 

^■^-0.4 


9 c+0.6 
^•^-0.5 


AiJ. I iS-C V f 




n 91 -1- n (14 


ns+°-"^ 

U.UO_Q 02 


R°° (km) 






IS*" 


Distance (kpc) 




8.4'' 


9 4-1-5.6 


TpSR 




2.1 ±0.8 


2.1 ± 0.8 


PpwN 




l.Sl'^ 


1.46'" 


/bt/nsa (10-" erg cm 


-2s-l) 


1.2 ±0.3 


i-4t^:^3 




-2s-l) 






/pj'" (10-"ergcm-2 
junabs (10-14 erg cm- 


s-1) 


2.3 ±0.4 


2.3 ±0.4 


-2s-l) 


^+6 
^-2 


5+5 
3-2 


/pwN (10-"ergcm- 
/PWN (10-1* erg cm- 


2s-l) 


2.8'" 


2.7*" 


-2s-l) 


4.0'" 


4.3*" 


Xi^/dof 




0.97/161 


0.97/161 



Note. — All uncertainties are 90% confidence inter- 
vals (i.e., 1.6(t). and /""''''= are the absorbed and 
unabsorbed fiuxes, respectively, in the 0.5-8 keV energy 
range. 



- A B-field strength of lO^^ G is assumed in the NSA 
model. 

*" - Held fixed during the fit. 

but with a significantly lower surface temperature of 
0.08lg'Q2 keV. The flux normalization suggests a distance 
of 2.4l^ g kpc, which is lower than the source distance of 
8.4 kpc. Table [2] summarizes the best-fit parameters of 
the BB+PL and NSA+PL fits. Afl uncertainties quoted 
are at the 90% confidence level. T he parameters are fully 
consistent with previous studies (jGonzalez et al.l 120051 : 
ISafi-Harb fc Kumail l200l) . and are better constrained 
in our case. For the case of magnetars, BB+ PL spec- 
tra are often observed (e.g., iNg et all 120 lit ) and the 
PL component is interpreted as magnetospheric upscat- 
tering of therma.1 pho tons from the surface emission 
(jThompson et al.]l2002t ). To check if this is the case for 
PSR J1119— 6127, we fitted the pulsar s pectra with a res - 
onant cyclotron scattering model (RCS: lRea et al.|[2008D . 
but obtained very poor results (x^ « 2), mainly because 
the model provides too little fiux above 3 keV to account 
for the observed hard X-ray tail. 

3.3.2. Phase-resolved Pulsar Spectroscopy 

To understand the nature of the X-ray pulsations, 
we performed phase-resolved spectroscopy using the PN 
data. Based on the pulse profile in Figure [2l we extracted 
the off-pulse spectra from phase 0.25-0.75, and the on- 
pulse spectra from the rest of the phase range. There 
are 673 ± 29 and 385 ± 23 net counts for the on- and 
off-pulse emission, respectively, in the 0.5-8 keV energy 
range. The spectra were grouped with at least 20 counts 
per bin, and fitted to the BB+PL model in the 0.5- 
8keV range. As shown in Figure [HJa), the pulsar emis- 
sion above 2.5 keV is dominated by non-thermal emis- 
sion, which shows no pulsations (see Figure [2|). There- 
fore, the non-thermal components (from both the pulsar 
and the PWN) were held fixed at the phase-averaged 
values during the fits. We also fixed the absorption col- 
umn density at the best-fit phase-averaged value of A'h = 
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Energy (keV) 

Fig. 3. — Top: best-fit blackbody plus power- law model for the 
phase-averaged spectrum of PSR J1119— 6127. The thermal com- 
ponent is shown by the dotted green line and the non-thermal com- 
ponent, which composes of pulsar and PWN emission, is shown by 
the dashed blue line. Only the FN spectra are shown for clarity. 
Bottom; on-pulse (red open circle) and off-pulse (black filled circle) 
FN spectra of the p ulsar, fitted with the blackbody plus power-law 
model. See Section 13.3.21 for details. 

2.2 X 10^^ cm~^. We note that the non-thermal flux is not 
enough to account for the off-pulse emission, suggesting 
the presence of thermal emission at pulse minimum. We 
found a higher BB temperature of kT = 0.22 ± 0.02 keV 
with an effective radius of i? = S.O^J gkm for the on- 
pulse emission. The off-pulse emission shows a hint of a 



lower temperature of kT — 0.15_q q3 keV and a larger ra- 
dius of 7^1^ km. The best-fit spectral models are shown 
in Figure [31[b) . We have also tried fitting the difference 
spectrum between on- and off-pulse emission. While it is 
consistent with a BB model, there are large uncertainties 
in the spectral parameters due to poor data quality. 

3.3.3. SNR Spectroscopy 

We extracted spectra of SNR G292.2-0.5 from the 
three brightest regions (see Figure [T]) from the XMM- 
Newton and Chandra data. Our region 1 roughly cor- 
respon ds to the "east" region in iGonzalez fc Safi-Harbl 
((2005I) . and regions 2 and 3 correspond to their "west" 
region. Point sources inside these regions were identified 
and excluded based on the exposure-corrected Chandra 
and XMM-Newton images. We note that region 2 is not 



covered by any Chandra observations, therefore, only 
the XMM-Newton spectra are used in this case. Sim- 
ilarly, only one Chandra exposure is useful for region 
3. Since the SNR covers a large part of the Chandra 
and XMM-Newton PN field of view, we estimated the 
background spectra using the blank-sky data p rovided 
by the calibration teams (jCarter fc ReadI I2007I ). which 
were taken from nearby regions in the Galactic plane. We 
also tried extracting the background spectra from source- 
free regions far off-axis in our observations and obtained 
similar results, although this way the sky background 
tends to be underestimated due to vignetting. After 
background subtraction, we obtained 1.1 x 10"* XMM- 
Newton counts and 5300 Chandra counts for region 1, 
4900 XMM-Newton counts for region 2, and 7900 XMM- 
Newton counts and 2000 Chandra counts for region 3 in 
0.5-7.5 keV range. 

Emission lines are found in all the spectra, albeit rather 
weak, indicating a thermal origin for the emission. For 
each region, we performed a joint fit to the spectra with 
all parameters linked, and we introduced a constant 
scaling factor between the XMM-Newton and Chandra 
spectra to account for the cross-calibration uncertainty 
(jTsuiimoto et al.l l 201lD. We first tried the s pectra l pa- 



rameters reported by IGonzalez fc Safi-Harbl |2005D , but 
these gave very poor fits. Using a non-equilibrium ion- 
ization thermal model fPSHOCK: lBorkowski et al.ll2001l) 
with solar abundances, we obtained reasonable results 
(xf, =1.3-1.5). The fits can be further improved by 
allowing the abundances to vary (with the VPSHOCK 
model) or by adding a hard component, and the im- 
provements are statistically significant as indicated by 
i^-tcsts. However, since the emission lines are weak, the 
VPSHOCK fits require rather low abundances (0.1-0.5) 
for all elements between Ne and Fe. Therefore, we believe 
that a two-component model could be a more physical de- 
scription of the thermal emission. In this case, the hard 
component can be described by either a non-thermal PL 
with photon indices F = 2.9-3.5 or by a high temper- 
ature {kT = 1-1.8 keV) PSHOCK model with very low 
ionization timescales (r < lO^scm"^) and solar abun- 
dances, and they provide the best fits overall {xl ~1.2). 
The results are listed in Table |3] and the best-fit two- 
temperature models are plotted in Figure SI In all the 
fits, the soft component has a large ionization timescales 
(r > 10^2 scm-^). The Chandra flux is ~ 30% higher 
than the XMM-Newto n flux, which is not atypical (see 
iTsujimoto et al.ll2011h . 

4. DISCUSSION 

In the following discussion, we adopt a source distance 
of 8.4 kpc and an age of 1.9 kyr. Although the true age 
of the system could be different than the pulsar's charac- 
teristic age, we stress that our discussion should remain 
qualitatively valid for a system age younger than a few 
thousand years. 

4.1. PSR Jl 119-6127 
4.1.1. Thermal Emission 

PSR J1119— 6127 is among the youngest radio pulsars 
with thermal emission detected. As indicated in Fig- 
ure [SJa), the pulsar flux below 2.5 keV is mostly ther- 
mal, and it is highly pulsed and shows good phase align- 
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TABLE 3 

Best-fit PSHOCK Plus Power-law (PS+PL) and Two-temperature PSHOCK (PS+PS) Models 

FOR THE Spectra of SNR G292.2-0.5 



Region 1 



Region 2 



Region 3 



Parameter 


PS-fPL 


py+py 


Pb+PL 


Py+Pb 


Pb+PL 


Pb+Pb 


A^H (10^2 


r, -1 7 + 0.07 
^- J- ' —0.08 


r, -1 7 + 0.10 

^•1 ' _0 08 


1 Q+0.5 


0.9 ±0.2 


1 ici+0.07 


1.0 ±0.1 


kTi (keV) 


0.18 ±0.02 


0.18+131 


^ — 0.02 


o.74i?,:JJ 


0.66«;«^ 


0.66 ± 0.03 


Ti (10l°scm~^) 




qn+90 


> 150 


> 30 


> 200 


> 100 


EMi (10^6 cm-3) 


360l?o° 


oor,+ 160 
ciOU_j20 


iinn+3000 


0.5 ±0.2 


1 9+0.4 


1 c+o.e 

-'•3-0.4 


r 


2.9 ±0.1 
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ment with the radio pulse (Figure [2]). Combining these 
two results suggests that the thermal emission could 
be dominated by a single hot spot near the magnetic 
pole. Further support for this picture is provided by 
the phase-resolved spectroscopy, which suggests a hint 
of a higher temperature for the on-pulse emission than 
the off-pulse emission (Figure M.^)). The BB+PL fit 
to the phase-averaged spectrum gives a high tempera- 
ture of 0.21 ± 0.04keV (= 2.4 ± 0.5 MK) with an emis- 
sion radius of 3^^ km, implying a bolometric luminosity 

of 1.9+0 g ^ lO'^'^ergs^^ for the thermal emission. The 
BB temperature and luminosity are among the highest 
of radio pulsars, even when compared to objects with 
similar age (see lOlausen et aT1[2010l: INg fc Kaspil 12011 
IZhuet al.ll20Tll) . Since the stellar surface temperature 
is non-uniform, the NSA fits only give the average tem- 
perature and it cannot be directly compared to the BB 
temperature. Indeed, the bolometric luminosity above is 
equivalent to that of a 13 km radius star with a uniform 
temperature of 0.096 keV (=1.1 MK), comparable to the 

b est-fit NSA temperatu re of 0.08l^-^| keV. 

IGonzalez et all (l2005l) and ISafi-Harb fc Kumaii (pOOSh 
pointed out that the BB luminosity and radius are too 
large to be reconciled with polar cap heating by return 
currents from the magnetosphere. Hence, the thermal 
emission is likely from cooling and its unusual properties, 
e.g. high temperature and large modulation, could be re- 
lated to th e stron g magnetic field of PSR J1119-6127. 
I Pons et all (j2007[) noticed an apparent correlation be- 
tween the BB temperature T and the dipole _B-field 
Bd of neutron stars, with T oc ^/Bd spanning over 
three orders of magnitude in Bd, from radio pulsars 
to magnetars. They attributed this to crustal heat- 
ing by magnetic field decay, and subsequently stud- 
ied the "magneto-ther m al ev olution" of neutron stars 
(|Aguilera et alJ l2008al lbl: |Po"ns et al.l l2009l ). In their 



model, heating from field decay increases the magnetic 
diffusivity and thermal conductivity of the crust, which 
in turn accelerates the field decay. As a result, neutron 
stars born with initial i?-fields > 5 x 10^^ G could have 
a surface temperature well above the minimal cooU ng 
scenario (|Aguilera et al.ll2008at [Kaminker et al.ll2009f ). 

The thermal emission of PSR J1119— 6127 offers an in- 
teresting test case for the above theory. While the BB 
temperature we o btained is comparable to the predicted 
value at the pole (jAguilera et al.|[2008aD . the bolometric 
luminosity of ~ 10"^^ — lO'^* erg s~^ could be eas ily ex- 
plained by passive cooling (e.g. iPage et all 120041) . thus, 
providing no evidence for energy injection from B-field 
decay. Indeed, the higher temperature at the pole could 
be attributed to anisotropic heat conduction. Since free 
electrons, which are the main heat carriers, are confined 
to the magnetic field lines, heat tr a.nsport perpendicular 
to th e field is strongly suppressed (|Greenstein fc Hartkd 
[1981 . This acts as a "heat blanket" at the magnetic 
equator, while the magnetic poles are in thermal equi- 
librium with the core. As a result, the surface tempera- 
ture distribution is non-uniform with warmer spots at the 
magnetic poles sep arated by cooler regions at lower mag- 
netic latitude (e.g. iGeppert et al.|[20M: iShabaltas fc Lail 
I20TI) . 

It is believed that the internal iTiagnetic field of 
a neutron star likely consists of a toroidal compo- 
nent, which could be generated by differential rotation 
that wraps the poloidal field in a proto-neutron star 
(e.g. [Thompson fc Dunca"^IT99l . IGeppert et all (poM l 
showed that the presence of a toroidal ZJ-field in the 
crust could further suppress the heat transport from the 
core towards the magnetar equator, because that requires 
crossing the toroidal field lines. As a result, this would 
enhance the heat blanket effect. Moreover, since the 
toroidal field component is symmetric about the mag- 
netic equator but the poloidal component is antisymmet- 
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Fig. 4. — Best-fit two-temperature PSHOCK models to the spec- 
tra of SNR G292.2— 0.5. The high- and low-temperature compo- 
nents are indicated by the dashed blue lines and dotted green lines, 
respectively. For clarity, only the PN spectra are shown and all 
spectra are grouped with a minimum of 50 counts per bin. 

ric, the total B-field would be asymmetric. Hence, the 
thermal emission from one magnetic pole could dominate 
over another, which could help explain the single-peaked 
pulse profile we found. Furthermore, the PF could be 
boosted by effects of limb- darkening and m agnetic beam- 
ing in the neutron star atmosphere (e.g., 'Pavlov et al.' 
[r994; .Raiagopal ct al..,1997; .van Adelsberg fc Lai.2006, ). 



Magnetic beaming is due to anisotropic scattering and 
absorption of photons in magnetized plasmas. The re- 
sulting radiation from a hot spot consists of a narrow 
pencil beam < 5° along the magnetic field and a broad 
fan beam at intermediate angles of 20°-60° (jZavlin et al.l 
1995; van Adels berg fc Lai 2006.) . 

To test these ideas, we employed a ma gnetized partiall y 
ionized hydrogen atmosphere model (|Ho et ahl |2008() , 
which accounts for the magnetic and relativistic effects, 
to generate X-ray spectra and pulse profiles with a sim- 
ilar procedure as [Ho (2007). We took a surface i3-field 
of 4 X IQ-^'^G and surface redshift of 1 + 2 = 1.19, and 
adopted the viewing geometry of magnetic inclination 
angle a w 125° and spin-axis inclination angle ^ « 145° 
obtained from modeling of the pulsar 7-ray lightcurve 
and radio polarization profile (Parent et al. 2011). Vac- 
uum polarization is ignored here, since it is not impor- 
tant for the total surface em ission at this field strength 
(jvan Adelsberg fc Lail 120061) . We tried different model 
parameters and found that two identical antipodal hot 
spots give a very low PF of ~ 5%. The observed single- 
peaked pulse profile can be produced if the temperature 
of one hot spot is higher than the other. Our best-fit 
model is given by a hot spot with effective temperature 
(observed at infinity) of kT = O.lSkeV (i.e. 1.5 MK) and 
an area of one third of the stellar surface, and for the 
other hot spot, we place a temperature upper limit of 
0.08 keV assuming a similar area. This is shown by the 
solid line in Figure [2j This temperature is lower than the 
best-fit BB temperature for the on-pulse spectrum (see 
Figure [3fb)) and also the emitting area here is larger. 
These are typical characteristics of an atmosphere model 
when compared to a simple BB, since non-gray opaci- 
ties in the former cause higher energy photons to be seen 
from deeper and hence hotter layers in the atmosphere 
(e.g. lRoman'illl987HShibanov et al.lll992D . Note that the 
model suggests a rather large emitting area. This is be- 
cause given the viewing geometry, the observer's line of 
sight would cut through the fan beam twice during one 
rotation. Therefore, a much smaller hot spot would pro- 
duce a double-peaked pulse profile. Finally, our model 
predicts a decreasing "max-min PF" with energy, from 
0.5 at 0.5 keV to 0.4 at 2keV, however, this difference is 
too small to be detected given our data quality. Deeper 
observations are needed to confirm this. 

As a caveat, we note that the viewing geometry ob- 
tained from the 7-ray lightcurve is model dependent. If 
the constraints on a and ^ are relaxed, the high PF could 
also be explained by some extreme view ing geometries, 
such that only one pole is observable (see lGonzalez et ahl 
l2005h . 

4.1.2. Non-thermal Emission 

The X-ray spectrum of PSR J1119-6127 clearly shows 
a PL component with a photon index of F = 2.1 ± 0.8 
and an unabsorbed flux of 5^2 x 10"^'* ergcm^^ s^^ in 
the 0.5-8 keV range. Since the RCS model gives a poor 
fit, this seems unlikely to be up-scattering of thermal 
photons as in magnetars. We believe that it could be 
synchrotron radiation originated from the PWN or from 
the pulsar magnetosphere. Assuming constant surface 
brightness, the PWN would provide an unabsorbed flux 
of only 8 x 10~^^ ergcm"^ s~^ in the 0.5-8 keV range 
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TABLE 4 

HiGH-B Radio Pulsars with Detected Pulsed Thermal Emission 



Name 


P 


B 


Tc 


Distance 


kT 


-^^bol 


PF 


Reference 




(s) 


(1013 G) 


(kyr) 


(kpc) 


(keV) 


(10^2 ergs^i) 






J1119-6127 


0.41 


4.1 


1.7 


8.4 


0.21 ±0.04 




0.48 ±0.12 


This work 


J1819-1458 


4.3 


5.0 


117 


~3.6 


0.12 ± 0.02 


1 2+2 


0.37 ±0.05 


1 


J1718-3718 


3.4 


7.4 


34 


~4.5 


0.19 ± 0.02 


4+5 


0.60 ±0.13 


2 



References. — (D IRea et aTll2009l : (2) IZhu et aL|[20Tll . 

Note. — The characteristic ages (tc) are inferred from the spin parameters. kT are the best- 
fit blackbody temperatures. The pulsed fractions (PF) are estimated with the "max-min PF" 
(see the text) in 0.5-2.5 keV, 0.3-5 keV, and 0.8-2 keV for PSRs J1119-6127, J1819-1458, and 
J1718— 3718, respectively. 



in the central 2'.'5-radius region. However, the actual 
PWN contribution is likely to be higher, since compact 
PWN emis sion is often found t o be peaked toward the 
pulsar (see iNg fc Romanl I2008D . Deeper Chandra ob- 
servations are needed to determine if the PL photon in- 
dices of the pulsar and PWN components are consistent. 
In addition to PWN emission, young rotation-powered 
pulsars generall y show strong m agnetospheric emission 
in X-rays (e.g. iNg et alll2008al) . Ho wever, this emis- 
sion is expected to be highly pulsed (jZhang fc Chend 
I2002f ). which is not observed in our case. At a distance 
of 8.4 kpc, the non-thermal emission we found has a lu- 
minosity of Lx — 4.2 X lO^^ergs"^ (0.5-8 keV), a few 
times larger than the flux 8.5 x 10'^^ ergs~^ suggested by 
iZhang fc Jiangl (j2005| ) based on the outergap model. 

4.1.3. Connection with Other High-B Pulsars and 
Magnetars 

Am ong all high-S radio pulsars listed in INg fc Kaspil 
(poll , only four (PSRs J1119-6127, J1819-1458, 
J1718-3718, and J1846-0258) are bright enough to have 
X-ray pulsations detected. Except PS R J1846-0258 
which has a purely non-thermal spectrum (jGotthelf et al.l 
dbOO; Ng et al. 2008a), the rest show thermal emission 
with high BB temperature around 0.1-0.2 keV. We sum- 
marize their timing and spectral properties in Table |4l 
Interestingly, their X-ray pulse profiles are all single- 
peaked and well aligned with the radio pulse, and show 
large modulations with PFs from 40% to 60%. These 
striking similarities suggest that the physics we discussed 
above in Section [4.1.11 could be applicable to other high- 
B pulsars. More X-ray observations are needed to con- 
firm this idea. In particular, PSRs J1734— 3333 and 
B1916-I-14, which are detected in X-rays but for which 
no pul sations have yet been found (^ ausen "eFaI1[20T2l: 
IZhu et al. 2009; Olauscn et al. 2010), would be ideal tar- 
gets. 

PSR J1846— 0258 is a very special object in the 
class, as it has shown m agnetar-like activity , includ- 
ing short X-ray bursts (|Gavriil et al.l 120081 ) accom- 
panied by a glitc h with unusual recovery properties 
(jLivingstone et alJ 120101 120 lit ). This led to specu- 
lation that high-B radio pulsars could be magnetars 
in quiescence. It has been proposed that radio pul- 
sars and magnetars could be a unified class of ob- 
jects (|Kaspi fc McLaughlinI l2005t |P erna fc PonsI [Ml 
and high-_B pulsars could exhibit magnetar-like bursting 
behavior, albeit less frequently (~1 per century). Ra- 



dio timing of PSR J1119-6127 reveale d a glitch in 2007 
that induced abnormal radio emission (jWeltevrede et all 
[201 If ). It is unclear if the pulsar exhibited any X-ray 
variability at the same time, due to the lack of sensi- 
tive X-ray monitoring. Our XMM-Newton observations 
in 2011 are not useful, since the post-outburst fiux re- 
laxation timescale is lik ely shorter than a few years (e.g. 
[Livingstone et aL|[2011|) . This highlights the importance 
of X-ray monitoring or rapid post-glitch follow-up in the 
study of high-i? radio pulsars and magnetars. 

4.2. SNR G292.2-0.5 
4.2.1. SNR Environment and Evolution 

We have carried out spectral analysis of three bright- 
est regions in SNR G292.2— 0.5. The spectra are well- 
described by absorbed two-component models. While 
the best-fit absorption column densities of the pulsar 
and region 1 are consistent (~ 2 x lO'^^cm"^), the val- 
ues are lower for regions 2 and 3 (~ 1 x lO^^cm"^). 
These seem to suggest an east-west gradient of A^h across 
the field. This is more evident in the three-color im- 
age in Figure [TJ the western part of the SNR is red- 
der than the eastern part, indicating a softer spectrum 
which could be due to a lower absorpt ion. This trend 
is consistent with previous find i ngs ([Pivovaroff et ahl 
200l [Gonzalez fc S afi-Harb I200l [Gonzalez et al.[ [20051: 
Safi-Harb fc Kumaii i200S ) . and can be attributed to a 
dark molecular cloud DC 292.3—0.4 to the east of the 
pulsar as first proposed bv ' Pivovaroff et all (j2001[ ). 

The SNR emission consists of a soft thermal compo- 
nent with kT < 1 keV and a hard component that can 
be described either by a non-thermal PL of photon index 
r = 2.9—3.5 or by high temperature non-equilibrium ion- 
ization plasma with kT =1 — 2 keV. While non-thermal 
synchrotron emission is often observed in young SNRs, 
the PL photon index we found (F ^ 3.5) is higher than 
the typical values of 2-3 in other systems (see a review 
bv [Revnoldd [20081 ). In addition, the fit statistics are 
slightly worse than that of the two-temperature model. 
These seem to suggest that the latter may be a more 
plausible model. In this case, the solar abundances of 
the high-temperature plasma suggest a circumstellar or 
interstellar origin, hence, the reverse shock probably has 
not yet propagated through the ejecta. Together with the 
high temperature and small ionization timescales, these 
provide support for the young age of the system. The 
emission measure could give us a handle on the ambient 
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density no- Assuming a slab geometry for the emitting 
regions with thickness equal to the width, and taking 
an electron-to-proton density ratio of 1.2 and a shock 
compression ratio of 4, we obtained hq ^ 0.1 cm~^ in 
all three regions. The low-temperature component, how- 
ever, is more difficult to interpret. Following the same 
calculation, we estimate rig of 0.6, 0.03, and 0.07 cm^'^ in 
regions 1,2, and 3, respectively. This could suggest dense 
clump for region 1. For regions 2 and 3, a comparison 
to the large ionization timescales (r > lO^^scm"'^) indi- 
cates that the plasma could have been shocked ~ 10^ yr 
ago, which seems difficult to reconcile with the SNR emis- 
sion. As we suggest below, the emission may arise from 
the progenitor wind interaction with the surrounding. 

At a distance of 8.4 kpc, the 8.'5-radius SNR shell has a 
physical size of 21 pc ( Caswell ct al. 2004). Such a large 
radius at a young age requires a high average expansion 
rate of Vs = 21pc/1.9kyr = 1.1 x lO'^kms"^. If this is 
the present-day shock velocity, then the proton temper- 
ature would be kTp = O.IItohWs = 130 keV, where mn 
is the proton mass, and an adiabatic index of 7 = 5/3 
and cosmic abundances (i.e., 90% H and 10% He) are as- 
sumed. A comparison to the best-fit kT « 1 keV implies 
an electron-to-proton temperature ratio of ^ 0.01, lower 
than in most SNRs observed (see Ivan Adelsberg et "aD 
|2008| ). If the ambient medium has a uniform density of 
no « 0.1 cm~^ as estimated above, then the swept-up 
mass would be ^ 100 Mq, much larger than the typical 
ejecta mass of < IOMq. Hence, the remnant could be in 
a transit ion to the adiab a tic ph ase. Based on the Sedov 
solution, iCrawford et al.l ()2001l ) derived a constraint on 
the ratio between hq (in cm~'^) and the explosion energy 
E51 in units o f 10^^ erg. The updat ed di stance and age es- 
timate s from lCaswell et aD (|2004l ) and I Weltevrede et al.l 
(|2011| ) imply E^i/uq = 200, requiring an exceptionally 
large explosion energy or a very low ambient density. 
It has been proposed that a newborn neutron star with 
magnetar-like i?-field (^ 10^^ G) and a short spin period 
{'^ 1 ms) would spin down quickly via magnetic brak- 
ing. As a result, the huge amount of rotational energy 
relea sed (^ 10^^ erg) could accelerate the SNR expan- 
sion C Allen fc Ho^^vath|[200l . While we cannot rule out 
this possibility, the explosion energies of three SNRs as- 
sociated with magnetars (Kes 73, CTB 109, and N49) 
are found to be close to the canonical value of 10^^ erg 
(jVink fc Kuipeiil2006[ ). which is not consistent with this 
hypothesis. 

As an alternative, the supernova could have occurred 
in a low-densit y w ind bubble (e.g., iGaensler et al.lll999l : 
iGvaramadze 200l iNg et al.ll2007D This scenario was 
briefly mentioned by iCaswell et al.l (l2004fl and here we 
further explore this idea. iChevalieii ()2005D discussed dif- 
ferent types of core collapse supernovae and the proper- 
ties of their remnants. Type IIP and IIL/b supernovae 
are the end products of red supergiants (RSGs), which 
have initial masses ^ 9 — 25 Mcr) and ~ 25 — 35M0, re- 
spectively (jHeger et al.l 120031 ). and they are surrounded 
by dense circumstellar RSG winds extending to a few 
pc. In contrast, supernovae of Type Ib/c are believed 
to originate from Wolf-Rayet (WR) progenitors, which 
are massive stars > 35Mq with high mass-loss rates and 
high- velocity winds, which arc about two orders of mag- 
nitudes faster than RSG winds. For a WR star evolved 



from an RSG, the fast WR wind would sweep up the slow 
RSG wind ejected earlier to evacuate a low-density wind 
bubble of radius > 10 pc ove r the lifetime ~ 2 x 10^ yr 
of the WR phase (iChevalied 120051 ) . The swept-up shefi 
is characterized by clumpy structure due to Rayleigh- 
Taylor instabilities, and ove rabundance in N and u nder- 
abundance in C and O (Gar cia- Segura et allll996[) . 

We argue that SNR G292.2-0.5 is unlikely to be a 
Type IIP or IIL/b supernova because in these cases 
the remnant should expand into the dense circumstel- 
lar RSG winds. A Type Ib/c event with a WR pro- 
genito r, one of the possibilities mentioned by IChevalieii 
(|2005l) . therefore seems more likely. In this picture, the 
SNR had a long free-expansion phase in the low-density 
wind bubble, with the initial expansion rate up to a few 
times 10^ km s~^, as seen in SN 1987A (Gae nsler et al.l 
ll997tlNg et al.l l2008bV When the shock reached the cav- 
ity boundary at large radius, it began to encounter the 
clumpy swept-up shell and hence decelerated. The shock 
interaction then gave rise to the observed hard compo- 
nent in the spectra, either thermal or non-thermal. This 
scenario may be able to explain the soft component too. 
The large ionization timescales of the emission are com- 
parable to the lifetime of the WR phase, suggesting that 
it may originate from the swept-up RSG wind. Two WR 
bubbles have been detected in X-rays and their spec- 
tra can be fitted by a t wo-temperat ure thermal model 
with fcT ~ 0.1 and 1 keV (IWri gge et al. 2005; Toala et al.l 
|2012| ). The high-temperature component is compara- 
ble to what we observed for SNR G292.2-0.5, while the 
low-temperature component may be too absorbed to de- 
tect in our case. In addition, hydrodynamic simulations 
and observa tions suggest a density of ~ 0.1 cm~'^ in the 
hot bubbles (Garcia-Se gura et al.lll996l : lDwarkadasll2007l : 
[^ala ct al. 2012), which is also consistent with our find- 
ings. Deeper X-ray observations in the future can probe 
the chemical composition of the circumstellar medium 
to confirm its nature. Moreover, if the SNR is detected 
in optical wavelengths, then spectral line observations 
could directly measure the shock velocity to reveal the 
evolutionary state of the remnant. 

4.2.2. Progenitors of High-B Pulsars and Magnetars 

In addition to PSR J1119-6127, only two other 
young high-i? radio pulsars have associated SNRs de- 
tected: PSR J1846-0258 in Kes 75 and PSR B1509-58 
in MSH 15—52. Intriguingly, b oth of them were 
claimed to have a W R progenitor ([Morton et al.l 120071 : 
IGaensler et al.l I1999D . There is also observational evi- 
dence hinting at a similar mass range for magnetar pro- 
genitors, e.g., IE 1048.1 -5937 CXOU J1647-4552, and 
SGR 1806-20 ( GacnslerelaU [20051 : IMuno et al.l [200fil : 
iFiger et al.l l2005f ). although intermediate masses were 
suggested for two cases: IE 1841-045 and SGR 1900-^14 
(IChevalied [2005t IDavies et al.l 120091) . This raises the 
important question regarding a possible correlation be- 
tween pr ogenitor mass an d neutron star magnetic field. 
Recently, lChevalie"rl ([201 1) compared different classes of 
young neutron stars in SNRs, including magnetars, ra- 
dio pulsars, and central compact objects, and suggested a 
tendency for higher i?-field objects to have more massive 
progenitors. However, it was noted that other param- 
eters, such as rotation, metallicity, and binarity, likely 
also play an important role. 
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The origin of magnetic fields in neutron stars has long 
been an open problem. Theories suggest that it could 
be the fossil B-fields of the progenitors preserve d dur- 
ing core collapse (e. g., Wolticr 1964; Rudcrnian 119721 : 
IFerrario fc Wickramasinghd |2006j) or from field ampli- 
ficati on via an a-fl dynamo in rapidly rotating stars 
(e.g., iDuncan fc Thompsonlll992t [Thompson fc DuncanI 
119931) . A positive correlation between progenitor mass 
and ne utron star .B-field may be expect ed in both sce- 
narios (jFerrario fc Wickrama singhe''2008"). Massive stars 
tend to have higher surface _B-ficlds (Petit ct al. 2008). 
They also spend less time in the hydrogen and helium 
burning phases, during which significant braking occurs. 
Hence, the core could retain a large angular momentum, 
resulting in efficient dynamo action. If further studies 
confirm that the progenitors of high-i? radio pulsars and 
magnetars have comparable mass, then no matter which 
of the above mechanisms is at work, the B-fields in these 
objects could be generated in a similar way. In this re- 
spect, magnetars could have the same formation chan- 
nel as radio pulsars, which would support a unificatioii 
of these classes of objects (Ka spi fc McLaughlinI [20051 : 
[Kaspi.,2010;,Perna k Pons 2011). 

5. CONCLUSION 

We performed a detailed X-ray study of the young 
high-S radio pulsar J1119— 6127 and its associated SNR 
G292.2-0.5 using deep XMM-Newton and Chandra ob- 
servations. The pulsar emission exhibits strong pulsa- 
tions below 2.5 keV, with a single-peaked profile that 
aligns with the radio pulse. Such a single-peaked pro- 
file and alignment seem common among thermally emit- 
ting high-i? pulsars, and we showed that the observed 
pulsed profile can be modeled by a large hot spot near 
the magnetic pole. The pulsar spectrum is well fitted by 
an absorbed BB plus PL model. The BB temperature 
kT = 0.21 ± 0.04 keV is highest among young radio pul- 
sars. The thermal emission has a bolometric luminosity 
of 1.9j^J g X 10'^^ ergs~^, which is consistent with neutron 
star cooling and no heating by B-field decay is needed. 
However, passive magnetic effects, including anisotropic 
heat conduction and beaming, must play a part in the 
high temperature and large modulation of the emission. 

The spectra of SNR G292.2-0.5 are best fitted by a 
two-component model with solar abundances. The soft 
component has a thermal origin with kT < 1 keV and 
a large ionization timescale, while the hard component 
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can be described either by a PL with F = 2.9 — 3.5 or by 
higher temperature {kT = 1 — 2 keV) thermal emission 
with low ionization. For a young age of 1900 yr inferred 
from the pulsar's spin down, the 21 pc SNR shell diame- 
ter (at a distance of 8.4 kpc) implies a fast expansion in 
the past. This could be the result of exceptionally large 
explosion energy > 10^^ erg or supernova in a wind cav- 
ity. We believe that the latter is more plausible, and this 
could suggest a Type I b/c supernova with a WR progen- 
itor. WR stars have a typical mass > 35 Mq and they 
have also been proposed as the progenitors for two other 
high-_B radio pulsars. There is evidence that magnetar 
progenitors could be in a similar mass range, leading to 
a speculation that the two classes of neutron stars may 
have the same formation channel. If confirmed, it will 
provide a strong support to the idea of unification of 
these classes of objects. 

We note that late in the pr eparation of this man uscript, 
we became aware of work ([Kumar et al.l l2012f) on the 
study of the SNR using the same archival Chandra data 
and a subset (^ 50 ks) of XMM-Newton data we have 
used in our study. While our deeper observations gave 
different spectral parameters for SNR G292.2— 0.5, we ar- 
rived at similar conclusions regarding a possible massive 
progenitor for PSR J1119-6127. 
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